Tube Width Fluctuations in F-Actin Solutions 
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We determine the statistics of the local tube width in F-actin solutions, beyond the usually 
reported mean value. Our experimental observations are explained by a segment fluid theory based 
on the binary collision approximation (BCA). In this systematic generalization of the standard 
mean-field approach effective polymer segments interact via a potential representing the topological 
constraints. The analytically predicted universal tube width distribution with a stretched tail is in 
good agreement with the data. 

PACS numbers: 61.25.H-; 82.35.Pq; 87.16.Ln 



The Edwards tube model provides a simple phe- 
nomenological description of the complicated topological 
constraints in entangled solutions of flexible polymers [1] . 
Using scaling arguments, Odijk [2 , Semenov [3] and oth- 
ers have generalized the idea to stiff polymers with a per- 
sistence length £ p larger than the characteristic backbone 
length between mutual collisions, which plays the role of 
the entanglement length in this context. Stiff polymers 
constitute an important and biologically relevant class 
of polymers, as they represent the major structural and 
stress-bearing units of the cytoskeleton of animal cells 
in the form of filamentous actin and microtubules [4]. 
Their large £ p and contour length L (on the order of 
10 |Ltm for actin) have opened the possibility of direct mi- 
croscopic visualizations of the tube [5]. Beyond its in- 
tuitive appeal and obvious relevance to single filament 
experiments in entangled solutions [TJ [8] , the tube enters 
explicit calculations of the overall mechanical properties 
of stiff polymer solutions [9HTT]. The latter have been 
measured rheometrically for biopolymers such as actin 
[T2HT5] and pectin [16] , as well as for self-assembling syn- 
thetic networks jTTJ [18] . While theoretical approaches so 
far have employed a homogeneity (or mean-field) assump- 
tion, treating the tube radius as a constant, experiments 
[SJ [7] EJ [19 , 20 indicate substantial heterogeneities, which 
recently also have been found in simulations [TO] . 

In this Letter, we present a systematic study of these 
heterogeneities in F-actin solutions, which we quantify in 
terms of the local tube radius profile R(s) (cf. Fig. [la- 
c). From R(s) we infer the entanglement length (FigTffl 
and, independently, the tube radius distribution P(R) 
(Fig. [I]i). The latter plays a similar role as the distribu- 
tion of void spaces or pores in other disordered materials, 
such as packings of grains or colloidal particles [21 . As a 
main result, P(R) is found to be described by a univer- 
sal master function with a stretched Gaussian tail (Fig.|3| 
inset). We develop a systematic theory of tube fluctua- 
tions that explains our data and provides the basis for 




c [mg/ml] 



FIG. 1: a) Superimposed confocal microscopy images of a flu- 
orescent actin filament in a background solution and a spline 
representing the tube backbone; scale bar: 5 |im; b) rectified 
image; c) tube radius profile R(s) determined as standard de- 
viation from Gaussian fits to the transverse intensity profile; 
d) normalized tube radius distribution P(R) obtained from 
cumulative contour lengths of 536, 804, 301, 225 and 116 (xm 
for actin concentrations c = 0.2, 0.4, 0.6, 0.8, 1.0 mg/ml, re- 
spectively. Solid lines represent a global fit by Eq. with 
p = 5.95 |xm _2 c[mg/ml] and L — 0.91 x L^°, including correc- 
tions accounting for the line spread function [6]. 



a more comprehensive characterization of the structural 
and elastic properties of stiff polymer solutions than pos- 
sible within the conventional mean-field theory [9]. 

F-actin solutions were prepared at various monomer 
concentrations c [6]. Rhodamine-phalloidin labeled solu- 
tions were mixed with unlabeled solutions at a ratio of 
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FIG. 2: a) A pair of colliding filaments is assigned a "distance" 
x and "topology" a — ± as exemplified for a = + with x > 
(left) and x < (right). A reflection at the horizontal mid- 
plane amounts to x — » — x, a — »• —a. b) Overlap region 
(shaded) for two segments of length L with orientations u, u' 
enclosing an arbitrary angle. 



1:1000. Time series of typically 150 pictures of individual 
labeled filaments were recorded using an inverse confocal 
microscope (LSM510, Carl Zeiss Jena, Germany). These 
pictures were superimposed to obtain a time-averaged 
image of the fluctuating test filament, where intensity re- 
flects the residence time of the filament (Fig.JlJi). Smooth 
tube contours connecting points of maximum intensity 
were constructed subject to a curvature- minimization 
constraint. The local tube radius R(s) was identified with 
the standard deviation of Gaussians fitted to the trans- 
verse intensity profiles of a rectified tube image, as ex- 
emplified in Fig. [T]o & c [6 . Along a single test fila- 
ment, R(s) exhibits pronounced fluctuations, which were 
binned to sample the tube radius distribution P(R). The 
result is depicted in Fig. [l]i (shaded areas) for various 
actin concentrations c. The peak position and width, 
corresponding to the typical value and the fluctuations 
of R, respectively, are seen to decrease with increasing c. 
Yet, as we demonstrate below, the non-Gaussian skewed 
and leptokurtic shape of P(R) is well described by a c- 
independent master function (Fig. [3j inset). 

We develop our theoretical approach along the lines of 
the BCA [9j [10] . It replaces the complicated topological 
constraints in an entangled polymer solution by an effec- 
tive model: an individual test polymer of length L in a 
tube-like harmonic confinement potential </>h 2 /2 per unit 
length, where h(s) parametrizes the transverse contour 
undulations along the backbone. In a self-consistent pair 
approximation, the tube potential is required to repre- 
sent the cumulative effect of independent pair collisions 
with the background polymers. These collisions are not 
to be understood as bare microscopic encounters of poly- 
mer backbones, though, but rather refer to effective tube 
collisions. The bare interactions caused by the impen- 
etrability of the polymer backbones are coarse-grained 
over the fast local conformational fluctuations of two col- 
liding polymers confined to their tubes, which prescribe 
a "quenched" geometry for them. The corresponding av- 
erages in a given tube configuration and over different 
tube configurations are represented by brackets (• • • )^ 
and over-bars 7TT , respectively. The over-bar on the tube 
stiffness <j) thus indicates that the latter represents the 



cumulative contribution from collisions in all possible ge- 
ometries and topologies as opposed to collisions in a pre- 
scribed tube configuration. A subtlety in counting states 
is that the topology of stiff but bendable polymers, as 
opposed to rigid rods, is not uniquely determined by the 
center-of-mass positions and orientations. Conversely, in 
any given topology a = ± (or "above" and "below"), 
positive and negative "distances" x ^ need to be dis- 
tinguished, as sketched for a = + in Fig. [2^i. 

The wormlike chain model in the weakly-bending rod 
limit with eigenmodes h(q) of the undulations h(s) is 
employed. Equipartition yields (h(q) • h(q f ))^ = 25 (q + 

q')/(Z P q 4 +?) for L > , where Lf = (£ p /^) 1/4 is the 
mean- field entanglement length [9]. Natural units with 
fc#T = 1, so that £ p is synonymous with the bending 
rigidity, are used throughout. The mean- field (projected) 
tube radius then follows as R°° = R(4>), given by the 
function 

r\4>) = |^<h 2 >, = 2- 3 /v /4 <r 3/4 (i) 

evaluated at <ft = (j). Similarly, one gets the parti- 
tion sum for the Gaussian contour undulations h(s) 
[9]. Its negative logarithm is the mentioned coarse- 
grained interaction or "BCA potential" F a (x) = 
-ln{erfc[-2- 1 / 2 ax/(^ + ^) 1 / 2 ]/2} between two tubes 
of radii Ro and R\ at separation x along the direction of 
nearest approach. If we allow for a uniform transverse 
displacement h of the test tube at an angle i(j relative to 
the ^-direction, which does not change the topology, the 
contribution to the confinement strength resulting from 
collisions in the prescribed tube configuration follows as 
cj) = L~ 1 d'^F a (x — hcosip)\h=o- In mean- field approxima- 
tion, setting Rq = Ri = R°° and performing the conflgu- 
rational average [9] yields the tube strength <p as a func- 
tion of R°° and the line concentration p (polymer length 
per volume). From Eq. 0, R°° = 0.66 x p- 3 ^£~ 1/5 is 
finally self-consistently determined [6]. 

The conventional BCA, as a mean-field theory, is ex- 
clusively concerned with the average values <\> and R°° . 
To get hold of the measured spatial tube width fluc- 
tuations, we introduce a canonical ensemble of N + 1 
independent entanglement segments of length L charac- 
terized by their individual fields <pi and corresponding 
tube radii Ri before averaging over the segment ensem- 
ble. In a formal generalization of the BCA that we call 
segment fluid model, any overlapping pair in the ensem- 
ble interacts with the BCA pair potential, now written 
as V(j(x) = x( r ? u ? u')F a (x) with the characteristic func- 
tion x of the overlap between two segments with orien- 
tations u, u' separated by r. As depicted in Fig. [2}d, 
the two segments are said to "collide" or to "overlap" if 
the projection of the center-of-mass of the segment with 
orientation u' onto the u-u' plane falls into the shaded 
area with edges of length L and the segment with ori- 
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entation u at its center. All pairs of segments are as- 
signed a binary topological state variable = ±. The 
confinement potential for a test segment with index is 
then computed as the cumulative effect from the colli- 
sions with all overlapping segments. Explicitly, after av- 
eraging over the uniformly distributed angles iJjoj of per- 
turbing displacements of the test segment, its individual 
tube strength is written as a local field 0o = J2iLi ^Oi 
with (j)Qi = (2L)~ 1 V" Q .(xoi). The average tube strength 
(j) is then obtained by taking the configurational average, 
i.e. by summing 0o, weighted by the Boltzmann factor for 
all pair interactions, over all topologies and positions 
and orientations of the segments. 

We now apply this formalism to calculate the distribu- 
tion P{(j>) of tube strengths. For the average over topolo- 
gies we exploit the identity e~ v+ + e~ v ~ = 1 for mutu- 
ally overlapping segments, such that Boltzmann factors 
not involving the test segment cancel out in the parti- 
tion sum. This is a consequence of the topological origin 
of the effective pair interactions. Moreover, due to the 
pair approximation, all interacting pairs are equivalent so 
that only one representative collision needs to be consid- 
ered explicitly. In the thermodynamic limit TV — »• oo, the 
characteristic function P(t) = e lt $ of P{(j>) then follows 
from the identity (1 + x/N) 
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P(t) = exp 



t<f>01—Vtr(xoi) 



(2) 



with the segment concentration n = (N + 1)/V = p/L. 
By a numerical inverse Fourier transform, P((f>) is ob- 
tained, which we recognize as the Holtzmark local field 
distribution of uncorrelated particles [22]. The tube ra- 
dius distribution P(R) then follows from P(4>) by ap- 
plying Eq. 0. In Fig. [3] (solid lines) it is shown for the 
special case of a pair potential V a with a common average 
range R4 = R for all segments i. 

To validate the result, numerical (Monte Carlo) inte- 
gration of a fluid of effective segments under the same 
assumptions as used in the theory was performed [6] , see 
Fig. [3] (dots). An examination of the numerical results 
establishes that the underlying asymmetric distribution 
P(<fi) is well approximated by a Gamma distribution 

r/e,6>( ( / ) )- Its shape and scale parameters, k = 4.013 x pLR 
2 

and 6 = 0.125 x (LR ) _1 , are determined by matching 
its first two cumulants with the predictions from Eq. ([2]). 
Via Eq. 0, an analytical expression for P(R) ensues, 



P(R) 



3R(k-l)\ 



exp(-y)y k , y 
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(3) 



It is compared with our numerical results for a fixed pre- 
scribed segment length L and various reduced concentra- 
tions pLR in Fig. [3] 

What remains to be done, is to identify the physi- 
cal meaning of the effective segment length L. Qual- 
itatively, one expects it to be equal to the mean-field 




FIG. 3: Dimensionless tube radius distribution P(R) (PdR = 

PdR) with scaling variable R = L~ 3/8 R~ 3/4 £l /8 R for re- 
duced concentrations pLR = 0.5,1,2 (bottom to top): 
Holtzmark distribution, Eq. (pi) (solid); Gamma distribu- 
tion approximation, Eq. Q (dashed); numerical integra- 
tion (dotted). A bimodal structure develops if the theory 
is pushed towards the (unphysical) limit pLR — > 0. In- 
set: Rescaled distribution as a function of R/R in semi- 
logarithmic representation: experimental data for c — 0.2 
(circles), 0.4 (squares), 0.6 (diamonds), 0.8 (upright trian- 
gles) and l.Omg/ml (downward- facing triangles), analytical 
approximation from Eq. (|3| as in Fig. [I] (dashed line), and self- 
consistent theory described in the main text with a slightly 
renormalized value of the segment length L — 1.62 x L e (solid 
line) . 



entanglement length = (£ p / '(j)) 1 ^ . The latter may 
in turn be expressed in terms of the mean-field tube ra- 
dius R°° and the line concentration p, respectively. Via 
Eq. 0, Lf = v / 2(^°°) 2/3 4 /3 oc p~ 2/5 . For the natural 
choice of L = the mean tube radius R = R((j)) fol- 
lows from a numerical solution of the implicit equation 
R JdRRP(R) a,sR= 1.15 x R°°. Its close match 
with the mean- field value R°° = lim^^oo R suggests 
to parametrize the distribution P(R) by the mean-field 
value R°° oc p -3 / 5 for i?, with negligible error. The shape 
parameter k thereupon becomes independent of p and 
only depends on L/L^° . The corresponding predictions of 
Eq. (|3| compare favorably with the measured tube radius 
distribution, as demonstrated in Fig. [l]i. In this compar- 
ison, L/L^° and p/c are used as global concentration- 
independent fit parameters. While the fit does indeed 
corroborate the expectation L/L^ w 1, p/c turns out 
to be about a factor of 6.6 smaller than estimated from 
the molecular weight and structure of monomeric actin 
[9]. This numerical discrepancy can be eliminated with- 
out significant consequences for the quality of the fit and 
the value oiL/Lf [13]. 

Note that Eq. Q implies that the concentration depen- 
dence enters the tube radius distribution only through 
the average tube radius i?, such that RP(R) defines a 
concentration-independent master function of R/R. The 
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FIG. 4: Left: Autocovariance function (ACF) of the tube 
radius profile averaged over 24, 37, 18, 13 and 9 filaments for 
c = 0.2, 0.4, 0.6, 0.8, and l.Omg/ml, respectively (symbols 
as in Fig. p] inset), with the abscissa rescaled by their initial 
slope A -1 (obtained by quadratic extrapolation for a varying 
fit interval s = so . . . s m ax, s m ax —> [6]). Because of the 
finite optical resolution of the microscope, data points below 
so = 230 nm were not considered. Right: A versus act in con- 
centration c with best power-law fit (exponent —0.44 ±0.09). 

inset of Fig. [3] (symbols) demonstrates that the data 
indeed scale satisfactorily. The semi-logarithmic repre- 
sentation reveals some systematic deviation of Eq. ([3| 
(dashed) from the data, however. This shortcoming 
is due to the preaveraging Ri = R employed in the 
derivation of Eq. and can be overcome by eval- 
uating our systematic theory more accurately. To 
this end, the variable tube radius Rq of the test seg- 
ment in the BCA potential is self-consistently identi- 
fied with the argument R of the distribution P(R), 
and only the radius R\ of the representative collision 
partner is preaveraged. This amounts to replacing R 
in Eq. Q by 2~ X I 2 \R 2 + (#°°) 2 ] 1/2 , before substitut- 
ing p = (R°° /0.66)~ 5 / 3 £ p 1 ^ 3 . Moreover, the segment 
length L is taken to be proportional to the local entan- 
glement length, L e = ^/2R 2 ^ 3 £p^ 3 . The resulting nor- 
malized tube radius distribution (solid line in the inset 
of Fig. [3| is in excellent agreement with our experimen- 
tal data. It exhibits (up to a spurious logarithm) a tail 
P(R) oc exp[— (R/R) 5 / 3 ], similar to what has recently 
been proposed on empirical grounds [20] . 

A further consistency check for the developed theory 
is provided by the spatial autocorrelation of the tube ra- 
dius profile in Fig.fTJ). It should decay over a characteris- 
tic length scale comparable to the entanglement length. 
The rescaled autocovariance functions determined from 
a large number of tube radius profiles for various concen- 
trations are shown in Fig. [4] (left panel), rescaled with 
their initial slope A -1 . Indeed, the concentration scaling 
of A thus obtained (right panel) compares favorably with 
the expectation A ~ ex c -2 / 5 . 

The fact that heterogeneities give rise to stretched tails 
in P(R) underscores the importance of shifting the atten- 
tion from characteristic "numbers" for R, L e etc. [12] to 
their skewed leptokurtic distributions. As demonstrated 
above, these are readily accessible in our BCA-based seg- 
ment fluid model. For instance, subtle non-steric correc- 
tions to the value of the tube radius R, e.g. due to tran- 
sient electrostatic attraction mediated by divalent coun- 



terions [24] , could systematically be studied by means of 
P(R) in the future. Our combined experimental and the- 
oretical results might also hold the key to a microscopic 
explanation of the ubiquitously observed broad distribu- 
tion of microrheological plateau moduli p~9j [25] . 
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